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ABSTRACT 

Using very high-resolution (R~ 125,000) and high quality (S/N > 350) spectra, we 
have searched for 6 Li in stars hosting extra-solar planets. From detailed profile-fitting 
of the Li I resonance line at 6707.7 A, we find no significant amount of 6 Li relative 
to the 7 Li for any of 8 planet bearing stars ( 6 Li/ 7 Li < 0.0 - 0.03 ) with a strong Li I 
lines. In particular, we do not confirm the presence of 6 Li with 6 Li/ Li = 0.13 reported 
by Israelian et al.(2001) for HD 82943, a star with two known planets. Several of the 
8 stars plus HD 219542 A, the planet-less primary of a binary, have been identified 
in the literature as possible recipients of accreted terrestrial material. For all of the 
planet-hosting stars and an additional 5 planet-less stars, we find no 6 Li. 
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1 INTRODUCTION 

Observers studying the abundance of lithium in stars con- 
tinue to uncover novel results. A recent novelty concerning 
main sequence stars was reported by Israelian et al. (2001) 
from measurements of the lithium isotopic ratio for HD 
82943, a star known to host giant planets, and HD 91889, a 
star similar to HD 82943 but lacking orbiting planets. For 
HD 91889, analysis of the 6707 A Li I resonance line showed a 
complete absence of 6 Li ( 6 Li/ 7 Li = -0.002 ± 0.006, by num- 
ber) in agreement with theoretical expectation that 6 Li but 
not 7 Li is thoroughly destroyed by protons in the pre-main 
sequence phase when the star possesses a deep convective 
envelope. The abundance of lithium, that is 7 Li, at log e(Li) 
= 2.52 is a typical value for a late F-type main sequence 
star. In striking contrast, HD 82943, a star with a similar 
total lithium abundance as HD 91889, showed a surprising 
amount of 6 Li: the ratio determined was 6 Li/ 7 Li = 0.126 
± 0.014. This extraordinary amount of 6 Li was cited by Is- 
raelian et al. as probable 'evidence for a planet (or planets) 
having been engulfed by the parent star'. Since lithium 6 Li 
and 7 Li - would presumably have been preserved in the un- 
fortunate planets, their ingestion added 6 Li (and 7 Li) to the 
atmosphere cleansed of 6 Li in the pre-main sequence phase. 

In this paper, we report measurements of the lithium 
isotopic abundance ratio for 8 stars (including HD 82943) 
known to host extra-solar planets, and for several similar 
stars not known to have extra-solar planets. Our principal 



result is that 6 Li is not detected in HD 82943 or other stars. 
The next section describing our observations is followed by 
a thorough description of our method of analysis. We use 
HD 82943 as one of the stars with which to illustrate our 
method. After presenting our results for the 6 Li/ 7 Li ratio, 
we discuss the implications of the non-detections of 6 Li, as 
well as offering a few remarks on the lithium abundance in 
the atmosphere of stars which host extra-solar planets. 



2 OBSERVATIONS 

The observed stars listed in Table 1 fall into two categories: 
stars with extra-solar planets, and similar stars for which 
there are no reported planets. Chosen stars have a Li I 6707 
A line of a strength enabling a meaningful limit (or detec- 
tion) to be set on the lithium isotopic ratio. Sharp-lined 
stars were given priority. The star 16 Cyg B, a host for 
planets but almost lacking in lithium, was observed in order 
to check on lines that might blend with the Li I 6707 A fea- 
ture. Comparison stars of similar atmospheric parameters 
were observed including HD 91889 with others chosen from 
surveys for lithium abundances by Balachandran (1990) and 
Chen et al. (2000). 

All of our observations were obtained with the coude 
echelle cross-dispersed spectrograph 2dcoude at the W.J. 
McDonald Observatory's 2.7m Harlan J. Smith reflector 
(Tull et al. 1995). The spectrograph was used in its higher 
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Table 1. Atmospheric parameters and observational details of the program stars. The absolute magnitudes, 
ages, and the masses were estimated using Hipparcos parallaxes and the apparent magnitudes combined 
with the isochrones and the evolutionary tracks computed by Girardi et al. (2000). 



Star 


Toff 


log 9 


6 


[Fe/H] 


M v 


Age 


M/Mq 


Rel|<j 




(K) 


(cm s z ) 


(km s _1 ) 






(Gyrs) 




HD8574 


6200 


4.20 


1.32 


0.12 


4.00±0.50^| 


2.0±2.0 6 


1.15 


PS 


HD 10697 


5605 


3.96 


0.95 


0.16 


3.73±0.05 


7.1±1.0 


1.10 


GG 


HD 52265 


6189 


4.40 


1.30 


0.28 


4.06±0.05 


1.0±0.5 


1.20 


GG 


HD 75289 


6140 


4.51 


1.47 


0.28 


4 n'S+n 04 


1.0±0.5 


1.20 


SIM 


HD 82943 


6010 


4.62 


1.08 


0.32 


4.35±0.05 


1.0±0.5 


1.00 


SIM 


HD 89744 


6338 


4.17 


1.55 


0.30 


2.79±0.06 


2.0±0.5 


1.40 


GG 


HD 141937 


6150 


4.30 


1.40 


0.31 


4.20±0.50 6 


0.8±2.0 f > 


1.22 


PS 


16 Cyg B 


5685 


4.26 


0.80 


0.07 


4.55±0.02 


9.0±1.0 


1.00 


LG 


HD 209458 


6063 


4.38 


1.02 


0.04 


4.29±0.10 


4.0±1.0 


1.10 


GG 


HD 219542 a| 


5989 


4.37 


1.20 


0.29 


4.50±0.50 6 


1.0±2.0 b 


1.15 


R.GG 


HD 75332 c 


6305 


4.49 


1.05 


0.24 


3.93±0.05 


l.Oil.O 


1.25 


GG 


HD91889 C 


6070 


4.41 


1.52 


-0.23 


3.75±0.03 


7.9±0.5 


0.95 


IS 


HD 142373 c 


5920 


4.27 


1.1 


-0.39 


3.62 ±0.02 


5.0±0.5 


0.90 


CYQ 


HD 154417 c 


5925 


4.30 


1.1 


-0.04 


4.46±0.04 


4.0±0.5 


1.05 


CYQ 


HD187691 C 


6088 


4.07 


1.35 


0.07 


3.66±0.03 


4.0±0.5 


1.25 


KF 



a References providing the adopted atmospheric parameters:- CYQ: Chen et al. (2000); GG: Gonzalez et al. 
(2001); IS: Israelian et al. (2001); KF: Fuhrmann (1998); LG: Laws & Gonzalez (2001); PS: present study; 
RGG: Gratton ct al. (2001); SIM: Santos et al. (2000). 

b Values of M t and Age are uncertain as stars lack either measured parallax or accurate visual magnitudes. 
c Stars not known to have planets. 



resolution mode. The recorded spectrum provides incom- 
plete coverage of 12 orders. The order containing 6707 A 
covers 25 A. The resolving power was R = A/AA = 125,000, 
as determined from the widths (FWHM) of thorium lines 
in the comparison spectrum. The spectra are over-sampled: 
there are 4.5 pixels per resolution element. 

Several exposures of a star were obtained with an ex- 
posure time of 15 to 20 minutes with total integration time 
of around 1 to 2 hours. A ThAr hollow cathode lamp was 
observed before and after a sequence of exposures on a star. 
Each image recorded on a Tektronix 2048 x 2048 CCD was 
processed in the usual way using IRAF subroutines to obtain 
a one dimensional spectrum. The wavelength calibration was 
applied to each individual exposure and was accomplished 
using a third-order polynomial fit to a large number of Th 
emission lines; this resulted in a typical rms accuracy of 2 
mA. 

Over the course of a sequence of exposures on a given 
star, the spectrum drifted on the CCD owing to the change 
in the Earth's motion, and to instrumental drifts. By cross- 
correlation, spectra were shifted so that they matched the 
first spectrum of a sequence. Continuum normalization of 
the resultant spectrum was achieved by fitting a sixth-order 
polynomial to each order. The S/N per pixel in the contin- 
uum near 6707 A is listed in Table 2. As a check on the 
continuum fitting, we compare the spectrum of HD 82943 
near 6707 A with that of the solar spectrum published by 
Hinkle et al. (2000) (Figure 1). It is seen that the line-free 
intervals in HD 82943 correspond very closely to those in 
the solar spectrum; the rms differences of about 0.2% are 
consistent with the S/N ratio of our spectrum of HD 82943. 



3 ANALYSIS: PRELIMINARIES 

In stellar spectra the Li I resonance doublet with a separa- 
tion of 0.15 A between the stronger blue and weaker red 
components is not resolved into two lines; the stellar line 
is an asymmetric blend with a greater width than other 
stellar lines of comparable strength. Addition of 6 Li with 
its stronger line almost at the wavelength of the weaker 7 Li 
line and its weaker line placed at an additional 0.15 A to the 
red introduces a wavelength shift of the blend and increases 
its asymmetry. Therefore, in order to extract the isotopic 
ratio from a stellar line, it is necessary to establish an ac- 
curate wavelength scale, and to quantify the various factors 
responsible for the line broadening. 

A necessary condition for defining the zero point of the 
wavelength scale is the determination of the stellar radial 
velocity. But examination of our spectra shows that it is 
not a sufficient condition because the radial velocity varies 
with the strength of the measured line. This is due to veloc- 
ity shifts induced by the convective motions (granulation) 
in the stellar atmosphere. The conventional list of factors 
influencing the width of a weak line like the 6707 A fea- 
ture includes the projected rotational velocity (vsini), the 
macroturbulence, the microturbulence, and the instrumen- 
tal profile. It is also necessary to consider the blends that 
contribute to the stellar 6707 A line. 

Reduced spectra were analyzed independently at the 
University of Texas (UT) and the University of Washington 
(UW). Results reported here for the full sample are those 
from the UT analyses, but are confirmed by the UW results. 
Some differences in the two approaches are noted later. 



0.9 



x 



0) 

> 

cti 

K 0.8 



0.7 




SUN 
• • • HD 82943 



6706 6708 

Wavelength (A) 



6710 



Figure 1. A comparison of the spectrum of HD 82943 and the Sun from 6705 A to 6710.25 A. Note the similar strengths and widths 
of the several Fel lines but the great difference in the Li I 6707.7 A feature. 



3.1 Radial Velocity and Convective Velocity Shifts 

The stellar radial velocity was measured for a sample of un- 
blended Fe I lines using the accurate laboratory wavelengths 
reported by Nave et al. (1994). There is a clear trend of 
velocity with equivalent width {W\) (Figure 2) in that the 
weaker lines are blue-shifted with respect to the stronger 
lines. This is a well known phenomenon attributed to gran- 
ulation - see the set of extensive measurements on Fe I lines 
in the Sun and Procyon reported by Allende Prieto & Garcia 
Lopez (1998) and Allende Prieto et al. (2002), respectively. 



The parallel with the Sun is almost exact, as might be ex- 
pected from Figure 1. Velocities of the same Fel lines have 
been measured from Hinkle et al.'s (2000) solar atlas (Figure 
2). 

Convective shifts are expected to depend on the atom or 
ion providing the line and on the line's excitation potential 
as well as its strength. In order to estimate the shift that 
might be expected of the Li I resonance line, we measured 
the solar shifts for lines in the red due to low excitation 
lines of Nal, All, and Cal: there are few of these lines in the 
limited bandpass of our stellar spectra. To a satisfactory 




precision, the velocity shifts of the additional lines follow 
the dependence of the Fel line shifts on W\. Since the Fel 
shifts for HD 82943 match well the solar shifts, we assume 
that the same equality holds for the Nal, All, and Cal lines 
in this star. It should be noted that expected differential 
velocity shifts between the Li I feature and metal lines in 
the same order are small (~ 0.2 km s _1 or less) relative to 
the isotopic shift of about 6.7 km s _1 but it is, as we note 
below, an observable effect. 

Granulation also renders stellar lines asymmetric. The 



line bisector is the usual way to express the asymmetry. 
Such an asymmetry is measurable for the stellar lines on 
our spectra. The bisector for a line in our program stars is 
similar to that of the same line in the solar spectrum. The 
asymmetry is very much smaller than the isotopic shift of 
the lithium line. 



3.2 Stellar Atmospheric Parameters 

Calculation of synthetic spectra requires a model atmo- 
sphere, and a line list. A model atmosphere was chosen 
from the grid provided by Kurucz (1995) according to the 
parameters T c g, log g, and [Fe/H] derived from published 
abundance analyses (see references given in Table 1. In this 
analysis we used an updated spectral analysis code MOOG 
(Sneden, 1973). The lithium isotopic analysis is insensitive 
to the particular choice of parameters. This and the fact 
that the published analyses are comprehensive mean that 
it is not necessary to derive afresh these parameters. For 
two stars with extra-solar planets, HD 8574 and HD 141937, 
with atmospheric parameters unavailable in the literature, 
the parameters T e g, log g, £ t , and [Fe/H] were derived using 
50 well defined Fel and 4 Fell lines with the required oscil- 
lator strengths ((//-values) for Fel lines taken from Nave et 
al. (1994) and from Giridhar & Ferro (1995) for Fell lines. 
We estimate uncertainties of about ±150 K in T e g, ±0.25 in 
log g, 0.25 km s" 1 in £ t , and ±0.10 in [Fe/H]. 

To match synthetic to observed spectra requires infor- 
mation on the various atmospheric effects that broaden the 
stellar profiles. The UT analysis considered three Fel lines 
in the same echelle order as the Li I feature and of a simi- 
lar strength to it. They are the lines at 6703.5, 6705.1, and 
6713.7 A. Synthetic profiles were computed using the appro- 
priate model atmosphere, the iron abundance, and the mi- 
croturbulence to match a Fel line's equivalent width. Two 
independent schemes were then employed to fit the synthetic 
line to the observed line profile. In the first scheme, the syn- 
thetic profile was broadened by a Gaussian modeled macro- 
turbulence defined by F 9 , representing the combined effects 
of stellar rotation, macroturbulence, and the instrumental 
profile as provided by the thorium lines in the comparison 
spectrum. In the second scheme, macroturbulence and rota- 
tional broadening were treated as different sources of broad- 
ening. 

In the first scheme, the T g was varied to match the 
observed profile of a selected Fe I line with the best fit judged 
by the reduced Xr where 

i 

where Oi and P» respectively are the observed and pre- 
dicted relative fluxes at data point i across the line profile, 
a is the rms error of the continuum, d = n — c (= 22 for 
Fel and 35 for Li I, typically) is the number of degrees of 
freedom in the fit, n is the number of observed data points 
involved in the fit, and c is the number of parameters to be 
determined from the computed line profile fit. Here, c is 3 for 
Fel lines (i.e., abundance, radial velocity, and the parameter 
r 9 ), and 4 for the Li I line with the isotopic 6 Li/ 7 Li ratio as 
the additional parameter. The run of \ 2 with T g is plotted 
and the minimum value found. 

The agreement between the predicted and the observed 
profiles is quantified with the x 2 -^est. The fit is considered 
very good when the \r — 1 f° r the best fit (Bevington & 
Robinson 1992). For HD 82943, Figure 3 shows the three 
Fel profiles and the run of Xr as a function of T g in the 
first scheme. Two of the lines are fit with a minimum Xr 
~ 1, and the third with a slightly higher value. The three 
results are consistent with a mean value F g — 5.03±0.12 km 



s 1 . Such a determination ignores the small line asymmetry. 
Rotational and instrumental broadening are subsumed into 

In the second scheme, all stars were also analyzed with 
vsini and V m as free parameters. The usini and V m values 
are determined for a given Fe abundance (determined from 
equivalent widths) such that the computed profiles best fit 
(i.e minimum Xi values) the observed Fe profiles. Figure 4 
shows the results for an Fel line for HD 82943 along with 
a similar star HD 75289. The three Fel lines were fitted in 
this way. The quality of the fits to the observed profile is 
quite similar to that obtained just using T 9 . Results for F g , 
Vm and vsini are given in Table 2 for all the program stars. 

In the UW analysis, the Fel lines at 5852.2, 6330.9, 
6703.5, and 6705.1, and 6715.4 A were employed to deter- 
mine the line broadening parameters vsini and V m . Instru- 
mental broadening was handled in the same manner as the 
UT group and the limb darkening coefficients were estimated 
using the work of Gray (1992). A linelist for 1 A regions sur- 
rounding each line was prepared via comparison of synthetic 
profiles generated using the MOOG spectral synthesis code 
(Sneden, 1973) at 0.001 A resolution with both the Solar 
Flux Atlas (Kurucz et al. 1984) and spectra of 16 Cyg B. 
As in similar prior work by Gonzalez (1998) and Takeda 
(1995), we discovered that in order to achieve good fits to 
these Fel features using accurate atmospheric parameters 
and accepted values for solar rotational and macroturbulent 
velocities, it was necessary to both reduce the microturbu- 
lent velocity, £ t , by a factor of 0.4 and allow for a 0.03 - 
0.05 dex offset in [Fe/H]. Applying this correction factor to 
£t results in a more accurate vsini. However, adopting a 
standard value of £ t would still result in a good fit to the 
line profiles. Synthetic profiles were calculated for each line 
using Kurucz (1995) model atmospheres, selected accord- 
ing to published values of T e s, log g, and [Fe/H] (Table 1); 
accompanying estimates of £t were adjusted in a manner 
consistent with that performed to fit the solar spectra. In 
each case, the precise placement of the continuum level was 
set to minimize the reduced Xi values of fits to adjacent 
continuum regions, and the velocity offset Av was adjusted 
to minimize the difference between the centroids of the ob- 
served and synthesized Fe lines. 

The UW analysis is based on 10-20 thousand synthe- 
sized profiles for each line at increments of 0.05 km s _1 in 
both vsini and V m and 0.01 dex in [Fe/H]. Values of Xi 
were calculated for each fit, resulting in 2-dimensional grids 
of Xr as a function of vsini and V m for 4-7 values of [Fe/H]. 
For each Fe I line, we adopt those values of v sin i and V m 
which yielded the minimum values of Xi amongst all com- 
peting fits. We take as our final solution for each star the 
average of the best-fit broadening parameters found from 
each line, and quote the standard deviation of their mean as 
an estimate of uncertainty. The final resulting broadening 
parameters for the UW group are given in Table 2. 

3.3 The Line List around 6707 A. 

Atomic data - wavelengths and gf- values - for the Li I 6707 A 
doublet and its isotopic and hyperfine splittings are of very 
high quality (Smith, Lambert, & Nissen 1998; Hobbs, Thor- 
burn, & Rebull 1999). For the UT line list, potential blend- 
ing atomic and molecular (CN) lines were identified from 
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Figure 3. Three Fe I lines in the spectrum of HD 82943 showing the fit of the synthetic profile for the optimum value of the Gaussian 
macroturbulence T 9 . The lower panels show the run of x? with T g . 



various sources: Kurucz (1995), Lambert, Smith, & Heath 
(1993), Brault & Miiller (1975), and Nave et al. (1994). The 
UT line list was tested and refined by matching the solar 
spectrum (Hinkle et al. 2000), and the spectrum of 16 Cyg 
B. Figure 5 shows the fit to the solar and 16 Cyg B spec- 
tra suggesting the line list lacks a line at about 6708.0 A. 
In their detailed analysis of the solar spectrum around 6707 
A, Brault & Miiller noted the presence of weak unidentified 
lines at 6708.025 A and 6708.110 A. Unfortunately, these 
lines straddle the weaker (red) line of the 6 Li doublet, and 
should, therefore, be included in the line list. On the basis of 



laboratory spectra available to them, Brault & Miiller con- 
cluded that these unidentified lines were not due to either 
CN or Fel. 

In earlier studies of the 6707 A line (e.g., Nissen et 
al. 1999), the presence of the unidentified lines was noted, 
but ignored as their contribution is negligible for metal- 
poor dwarfs. However, in F-G main sequence stars of ap- 
proximately solar metallicity, this pair of lines will signif- 
icantly affect the red wing of the 6707 A profile. King et 
al. (1997) noted the effect of these unidentified lines on the 
solar 6 Li/ 7 Li ratio. Treating these missing lines as 6 Li com- 
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Figure 4. Determination of v sin i and V m for HD 82943 and HD 75289 



ponents they obtain 6 Li/ 7 Li = 0.10 for the Sun. This high Li 
isotopic ratio for the sun is very unlikely for obvious reasons, 
and they attributed the lines to Si I. In the solar spectrum, 
the 6708.025 A and 6708.110 A lines have equivalent widths 
of W\ = 1.4 mA, and < 0.5 mA, respectively. Our analyses 
of the solar and 16 Cyg B spectra give lithium ( 7 Li) abun- 
dances consistent with previous analyses. For the Sun, our 
line list with the Kurucz solar model gives log e(Li) = 0.98 
and for 16 Cyg B we obtain log e(Li) = 0.71 which are close 
to the values of loge(Li) = 1.05±0.06 and loge(Li) <0.60, 
respectively, derived by King et al. (1997). 

Our analysis of selected stars was made first without 



the unidentified lines in the line list. Then, we repeated the 
analysis on the assumption that they were neutral Ti lines 
with a lower excitation potential of 1.88 eV and derived their 
<7/-values from a fit to the solar spectrum. As a check, we 
synthesised the spectrum of 16 Cyg B and found that the 
(unidentified) Til lines matched the observed feature well 
(Figure 5). The star 16 Cyg B is cooler than most of our 
program stars and is very deficient in lithium (King et al. 
1997). 

Figure 5 shows that while the 7 Li feature gets weaker, 
the missing feature at 6708.0 A is a little stronger in 16 Cyg 
B than in the solar spectrum. We have varied the assumed 



Table 2. Lithium isotopic 6 Li/ 7 Li ratio results derived from 6707 A Li-profile. One-tr errors are quoted. 



Star 


S/N 


L 9 

(km s -1 ) 


v sini 
(km s — 


J 


v m 

(km s — 


J 


v sinif 
(km s~ 


I) 

J 


m 

(km s — 


) 


logefLifl 

(dex) 




A 


2 b 
min 


HD 8574 


300 


7.22zt0.16 


3.30zt0. 


20 


4.65±0. 


10 


3. 


60zt0. 


.26 


5.30±0. 


.30 


2.71±0.05 


0.03±0. 


.05 


l 


.02 


HD 10697 


300 


5.27±0.10 


1.20±0. 


30 


4.30±0. 


15 


1. 


61±0. 


.19 


4.04±0. 


.18 


1.91±0.03 


O.OOztO. 


.06 


l. 


.15 


HD 52265 


450 


7.39zt0.08 


3.10±0. 


25 


5.50±0. 


.13 


3 


.88±0. 


.11 


4.81±0. 


.16 


2.76±0.01 


O.OlztO. 


.03 


l. 


.05 


HD 75289 


350 


7.08zt0.10 


3.40±0. 


25 


4.50±0. 


10 


4. 


,17±0. 


.25 


4.15±0. 


.30 


2.77±0.02 


0.02±0. 


.03 





.90 


HD 82943 


450 


5.03±0.12 


1.60±0. 


20 


3.65±0. 


10 


2. 


09zt0. 


.19 


3.73±0. 


.13 


2.43±0.02 


O.OOztO. 


.02 


l. 


.06 


HD 829430 


530 


5.00±0.14 


1.55±0. 


13 


3.73±0. 


15 












2.43±0.02 


O.OlztO. 


.03 


l. 


.20 


HD 89744 


550 


13.54±0.20 


7.60±0. 


28 


S.OOztO. 


20 


8. 


.02±0. 


.13 


7.71±0. 


.32 


2.11±0.01 


O.OOztO. 


.03 


l. 


.20 


HD 141937 


380 


5.14±0.25 


1.15±0. 


30 


3.85±0. 


15 


1 


89zt0. 


.26 


4.00±0. 


.29 


2.55±0.04 


O.OlztO. 


.03 


l. 


.12 


16 Cyg B 


600 


4.79±0.10 


1.60±0. 


25 


3.30±0. 


18 


1 


64±0. 


.16 


3.49±0. 


.16 


0.71±0.02 


O.OOztO. 


.03 


l. 


.30 


HD 209458 


150 


7.24±0.05 


3.33±0. 


15 


4.83±0. 


18 


3 


57zt0. 


.11 


5.24±0. 


.22 


2.70zt0.02 


O.OOztO. 


.03 


0. 


.96 


HD 219542 A 


300 


4.50±0.20 


1.55±0. 


25 


3.55±0. 


10 


2. 


.24±0. 


.16 


3.66±0. 


.12 


2.26zt0.03 


0.03zt0. 


.04 


1. 


.06 


HD 75332 


650 


12.62±0.21 


8.05±0. 


20 


6.78±0. 


23 


8. 


20zt0. 


.06 


6.60±0. 


.16 


3.18±0.01 


O.OOztO. 


.02 


1. 


.05 


HD 91889 


600 


6.63±0.12 


2.20±0. 


30 


5.15±0. 


16 


2. 


89zt0. 


.26 


4.49±0. 


.14 


2.46±0.02 


O.OOztO. 


.02 





.98 


HD 142373 


700 


5.12±0.13 


1.05±0. 


20 


4.28±0. 


15 


2 


.21±0. 


.51 


3.67±0. 


.46 


2.52±0.01 


O.OOztO. 


.01 


1. 


.36 


HD 154417 


600 


8.68±0.18 


3.95±0. 


32 


6.35±0. 


18 


4. 


90±0. 


.06 
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a Values independently derived by UW group. 

b Values computed using the single Gaussian broadening parameter r s . 
c Two spectra of HD 82943 were obtained and independently analysed. 



excitation potential of the Til lines, and tested to changes 
in Tea- We found a Til line with a LEP of 1.88 eV varies 
with Trff by around 0.08 dex/100 K, and for a LEP of 0.0 
eV, the Til line strength varies by 0.12 dex/100 K. Thus, 
by adopting 1.88 eV lines of neutral Ti for the unidentifed 
lines, instead of 0.0 eV lines, we overestimate their strength 
by around 0.10 dex for stars of T e g = 6400 K, the hottest 
stars in our sample. Higher LEP lines of Ti and lines of 
higher ionisation potential elements like Si and Mg also in- 
vestigated. Within our T e a range, the T e a effect on strength 
of these lines does not have noticeble effect on 6 Li/ 7 Li ra- 
tios. Given that our programme stars are quite similar to 
the Sun, identification of the pair as Ti I lines is not critical 
to the analysis, but their inclusion is critical. The final line 
list is given in Table 3. 



4 THE 6 Li/ 7 Li RATIO 

Addition of 6 Li to the mix of lithium isotopes shifts the 
6707 A line to the red and increases the line's asymmetry. 
For the meteoritic isotopic ratio ( 6 Li/ 7 Li ~ 0.10), the center- 
of-gravity of the profile shifts redward by ~ 710 m s _1 with 
respect to the 7 Li-only profile. Under ideal conditions, where 
the observed line shifts are independent of their strength and 
carrier, it is straightforward to determine the wavelength 
scale from Fel and other lines and to estimate the 8 Li contri- 
bution from the 6707 A line's wavelength. Convective shifts 
in the atmosphere introduce a trend in the relation between 
a line's wavelength and strength (Figure 2). In the case of 
HD 82943, the Fe I lines of strength W\ « 50 mA are shifted 
by 120 m s _1 relative to lines of strength W\ ~ 10 mA. The 
differential shift is evident from syntheses of the 6707 A re- 
gion. When all lines are given an identical radial velocity, 
it proved impossible to fit simultaneously the Li I line and 
the adjacent weaker Fel at 6707.433 A. While the latter line 



was matched in strength and width, the wavelength separa- 
tion between it and the Li I line was not. When the stellar 
radial velocity is measured from the Fe I lines at 6703 A and 
6705 A, a pair of lines similar in strength to the Li line, the 
observed wavelength of the deepest part of the Li I line is at 
the predicted wavelength for 6 Li/ 7 Li = 0.0 to within about 
1 mA. This predicted wavelength is moderately sensitive to 
small amounts of 6 Li, for example, the shift from 6 Li/ 7 Li = 
0.0 to 0.1 is 0.005 mA (~ 220 m s _1 ), but use of the wave- 
length is compromised by the presence of convective shifts. 

As our primary monitor of the lithium isotopic ratio, 
we use the profile of the 6707 A line. This is less affected 
by the granulation. The Li abundance and the 6 Li/ 7 Li ra- 
tios for the program stars are derived with the following 
procedure by the UT group. Using the model parameters 
given in Table 2 and the line list given in Table 3, Li-profiles 
are computed for a given 6 Li/ 7 Li ratio by varying the Li 
abundance and the wavelength of the Li I line such that the 
predicted profile fits best the observed profile. We also took 
the liberty to adjust the observed profile vertically within 
the rms error (~ S/N -1 ) of the continuum. Abundances of 
elements (C, N, Ti, V, and Fe) contributing lines that are 
blended with the Li-profile are taken from the sources which 
provided the atmospheric parameters (see Table 1 for refer- 
ences) . Predicted profiles are first broadened with the single 
Gaussian parameter F g and compared with the observed 
profile. Then, the procedure is repeated but the predicted 
profiles are broadened using the combination of vsini, V m , 
and the instrumental broadening. 

Sample results are shown in Figure 6 for HD 82943 for 
the case where rotation and macroturbulence are treated 
as separate effects. Synthetic profiles for the case of Gaus- 
sian broadening parameter r g alone are not discernibly dif- 
ferent. The agreement between the predicted and the ob- 
served profiles is quantified using Xr- The minimum values, 
Xmin> obtained from Xr analysis are given in Table 2. In Fig- 
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Figure 5. Observed and synthetic spectra around the Li I 6707 A feature for the Sun and 16 Cyg B. The contribution of the unidentified 
lines near the weaker 6 Li line is given by the difference between the dashed and solid lines between 6707.9 A and 6708.2 A. 



ure 6, the minimum \ 2 variation values, Ax 2 , for HD 82943 
obtained for different 6 Li/ 7 Li ratios are plotted. The A^ 2 
is computed by subtracting \i from the Xf ! where Xf is 
computed by dropping wavelength as a free parameter. The 
wavelength is fixed for the minimum Xmin value from the 
Xi analysis. Thus, by definition both the x? arid Xf should 
have the same minimum value of xLin- For HD 82943, we 
obtained the abundance log e(Li) = 2.435±0.005 and the 
ratio 6 Li/ 7 Li = 0.00±0.03. The quoted 1-a errors in the iso- 
topic ratios (Table 3) are obtained from the x 2 variation, 
the formal 1-a, 2-a, and 3-cr errors (Bevington & Robinson 



1992) are indicated in Figure 6. The quoted errors in the 
log e(Li) are 1-a errors obtained from the x 2 variation of 
the Li-profile fitting, and do not include errors arising from 
the choice of model atmosphere parameters, especially the 
effective temperature. 

In Table 2, we give the results for the lithium abundance 
and isotopic ratio obtained in the UT analysis using the as- 
sumed single line broadening parameter T g . Also, values of 
Xmin for the final best fit between the observed and the pre- 
dicted Li-profile are given in Table 2. The abundance and 
the isotopic ratio show very weak dependence on how the 




Figure 6. Computed synthetic profiles for different 6 Li/ 7 Li ratios are compared with the observed Li-profile. Computed profiles are 
broadened with the broadening parameter r 9 (top panel) and with the separated radial-tangential macroturbulcncc, rotational velocity 
and the instrumental profile (bottom panel). The \ 2 variance A\ 2 versus the Li isotopic ratio is shown for both the cases. The errors 
1-cr, 2-cr, and 3-cr are noted. 



velocity broadening is assigned. Additional examples of ob- 
served and synthetic spectra and accompanying plots of \ 2 
are shown in Figures 7 & 8. Although the lithium abundance 
is sensitive to the adopted model atmosphere - temperature 
is the controlling influence - the isotopic ratio is quite in- 
sensitive to it. We repeated the analysis using models drawn 
from the MARCS grid (Gustafsson et al. 1975) and found 
no noticeable difference in either the broadening parame- 
ters or the Li-isotopic ratio. Table 2 shows that, across the 



sample of stars, the lithium is pure 7 Li with no convincing 
detection of 6 Li. This result holds for stars with and without 
extra-solar planets. 

Our analysis of the 6707 A line profile using standard 
model atmospheres obviously neglects the influence of the 
stellar granulation, which is responsible for the line shifts 
and asymmetries. These can have only a slight effect on the 
results. Stellar equivalents of sunspots and faculae might 
cause the lithium line to be strengthened over a few localized 



Table 3. Adopted linelist at the vicinity of the Li 6707A profile 
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regions of a stellar disk. Consider, for example, a strength- 
ening of the line in an area close to the receding limb of the 
star. The result would be a red asymmetry to the lithium 
line, and quite possibly a weaker asymmetry would be im- 
posed on the Fe I comparison lines used to set the broadening 
parameters. In these circumstances, one might overestimate 
the contribution of 6 Li. Similary, if areas on the approach- 
ing limb carry undue weight in forming the profile, the 6 Li 
contribution could be underestimated. If a major portion of 
the Li I line is formed in a few special areas, our theoretical 
profiles would fail to fit the observed profiles. 

Our demonstration that excellent fits to the 6707 A pro- 
files are possible is a good indication that stellar spots and 
faculae cannot play an important role. Israelian et al. rec- 
ognized that a false 6 Li signal could arise from the neglect 
of the effects of localized regions on the formation of the 
6707 A line. They advanced cogent reasons - photometric 
and radial velocity stability - for rejecting this possibility in 
the case of HD 82943. The fact that our observed profiles 
are very similar to their profile (see below) is additional ev- 
idence that the lithium line is not formed preferentially in a 
few areas. 

4.1 The Case of HD 82943 

Israelian et al. (2001) derived the high value of 6 Li/ 7 Li = 
0.126 ± 0.014 for HD 82943. Our analysis (UT and UW) 
shows no indication of 6 Li or 6 Li/ 7 Li = 0.0±0.03. One can 
think of two possibilities for the discrepant results: a) a dif- 



ference in the observed spectra, and b) differences in the 
analyses. 

Israelian et al.'s spectrum is of similar quality to ours 
with respect to resolving power (R = 110,000 vs our 125,000) 
and S/N ratio in the continuum near the Li I line ( S/N~ 500 
in both the studies). Comparison of their electronically pub- 
lished Li I profile for HD 82943 with our profiles shows near 
perfect agreement; the differences are much smaller than 
those between synthetic profiles for 6 Li/ 7 Li = 0.0 and 0.13. 

Given that there is no difference between Israelian et 
al.'s and our spectra of a size sufficient to account for the 
different 6 Li/ 7 Li ratios, the sources of the discrepancy must 
lie in the methods of analysis. In seeking an explanation, 
we consider the adopted model atmospheres, the fitted line 
broadening parameters, and the chosen line lists. We adopt 
the model atmosphere parameters - T e g , \ogg, [Fe/H], and£ t 
- used by Israelian et al. They do not explicitly indicate the 
source of the model atmosphere but it is most improbable 
that given the similarities between modern grids of model 
atmospheres (see our test of a MARCS model) that an alter- 
native choice can lead to anything but a very small change 
in the recovered 6 Li/ 7 Li ratio. 

In obtaining the broadening parameters, Israelian et al. 
used two of our three Fei lines to obtain vsini = 1.65 ± 
0.05 km s _1 and V m = 3.90 ± 0.2 km s _1 . These com- 
pare well with our values vs'mi — 1.60 ± 0.20 km s~\ and 
V m = 3.65± 0.10 km s _1 , where we assumed the same form 
for V m as Israelian et al. Analysis of our spectrum using 
their broadening parameters also gives 6 Li/ 7 Li = 0.01, not 
0.126. Our result is insensitive to the form adopted (Gaus- 
sian macroturbulence vs rotation and macroturbulcnce) for 
the line broadening (Figure 9). 

The inferred contribution from 6 Li to the observed pro- 
file depends on the completeness and accuracy of the line 
list. It would appear from their paper that Israelian et al. 
did not include the unidentified line at 6708.025 A in their 
line list, and overlooked the weaker line at 6708.125 A. The 
unidentified lines which we assign to Ti I play a role because 
they lie on either side of the weaker 6 Lil line. Their com- 
bined strength in the solar spectrum is about 1.5 mA. In 
a metal-rich star like HD 82943, their strength is greater 
than in the Sun. Given that the equivalent width of the 
6707 A feature is 50 mA for HD 82943, the lines' estimated 
contribution of 2.2 mA to the weaker 6 Li line corresponds 
to 6 Li/ 7 Li ~ 0.06. The linelist compiled independently at 
the University of Washington differs very little from that in 
Table 3. Use of this linelist instead of that in Table 3 does 
not give significantly different synthetic spectra computed 
for the same input data. Israelian et al. used HD 91889 as a 
comparison star, but this is not necessarily a perfect foil for 
eliminating the effects of weak blends. This star is metal- 
poor with [Fe/H] = -0.23 or 0.55 dex down from HD 82943. 
This relative metal deficiency plus slight differences in effec- 
tive temperature and surface gravity mean that unidentified 
lines like the supposed Ti I lines will be appreciably weaker 
in HD 91889 than in HD 82943. In contrast, our sample of 
comparison stars includes examples of similar metallicity to 
HD 82943. 

To investigate further the influence of the unidentified 
lines, we analysed the HD 82943 profile with different Ti 
abundances assigned to the pair of unidentified lines. As the 
abundance is decreased, the drop in absorption from the Ti I 
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Figure 7. Determination of 6 Li/ 7 Li ratio by fitting the Li-profile for two stars with extra-solar planets: HD 75289 (top panel) and 
HD 52265 (bottom panel). The x 2 -test results are shown for both the stars. 



lines is compensated for an increasing Li abundance, but 
the quality of the fit to the observed line decreases steadily. 
This is well shown by Figure 10 where we plot the Ti abun- 
dance versus the derived 6 Li/ 7 Li that best fits the line. The 
exercise was done for Til lines of 0.0 eV and 1.88 eV. In 
the latter case, the x f° r the best fit is marked and clearly 
increases as the contribution of the unidentified lines is re- 
duced. 

Our stellar sample include several other metal-rich 
stars, including one not known to host planets. A plausible 
assumption is that, if 6 Li is present in these stellar atmo- 



spheres, the 6 Li/ 7 Li ratio will vary from star-to-star, and is 
likely absent in the star without planets (HD 75332) . Incor- 
rect representation of the list of blending lines could bias 
the 6 Li/ 7 Li determinations but is unlikely across this sam- 
ple of metal-rich stars to alter the spread in 6 Li/ 7 Li ratios. 
The fact that we find 6 Li/ 7 Li ~ 0.00 in all cases suggests 
that 6 Li is not present in these metal-rich stars, including 
HD 82943. Blends have a smaller effect on the more metal- 
poor stars. The best fits to the observed profiles are found 
for 6 Li/ 7 Li = 0.00 (i.e., non-negative) is a good indication 
that the blending lines have not been overestimated. 




Figure 8. Same as Figure 7. Li-profile fitting and the x 2 -tcst results for HD 219542 A (top panel) and HD 91889 (bottom panel) which 
are not known to have planets. The Li-profile of HD 219542 best fits the computed Li-profile of 6 Li/ 7 Li = 0.03 ±0.04. 



5 DISCUSSION 

Our sample of planet-hosting stars are placed in an HR- 
diagram (Figure 11) as a guide to their evolutionary state. 
We use the Hipparcos parallaxes, the observed V magnitudes 
(interstellar reddening is neglected), and the reported T e g. 
On this diagram, we place isochrones from Girardi et al. 
(2000) for [Fe/H] = 0.0 and 0.3 and ages from 1 to 9 Gyr. 
All but three of the stars are young; their ages are 4 Gyr or 
less (Table 1). 

The young stars have a lithium ( 7 Li) abundance that 



appears normal for stars of their age and metallicity. An em- 
pirical demonstration of normality is provided in Figure 12 
where the lithium abundance is plotted versus T e ff and com- 
pared with observations from three open clusters of differing 
ages - Pleiades at 0.1 Gyr, Hyades at 0.7 Gyr, and M67 at 4.5 
Gyr - summarised by Jones, Fischer, & Soderblom (1999). 
Ages inferred from Figure 11 are consistent with those listed 
in Table 1. Lithium abundances for our main sequence stars 
and the clusters are approximately consistent with theoreti- 
cal predictions (see, for example, ProfBtt & Michaud 1989). 
We note that our sample stars were selected to have mod- 




Figure 9. Li-profile fitting for HD 82943 using three different forms of line broadening parameters (left panel): F g = 5.03 km s -1 ), 
vsini = 1.60 km s _1 , V m = 3.65 km s _1 from the UT analysis and vsini = 2.09 km s -1 , V m = 3.73 km s — 1 from the UW analysis. 
In the right panel we showed x?-test analysis for each case. Note a very slight change in the quality of the fit but in all the three cases 
Xmin is found for 6 Li/ 7 Li < 0.01. 



erate to strong Li lines. So, the plot is necessarily biased 
towards higher Li abundances. 

Normal behaviour for 6 Li must be defined on theoreti- 
cal grounds because there are no observations of this isotope 
in young stars. Standard models (Profntt & Michaud 1989) 
predict effectively complete destruction of 6 Li in the pre- 
main sequence phase. Destruction results when the convec- 
tive envelope exposes lithium to warm protons at its base. 
Little to no destruction of 7 Li occurs; the (p, a) reactions 



occur about 80 times more slowly for 7 Li than for 6 Li. De- 
struction of 7 Li occurs on the main sequence, as revealed by 
the lowering of the lithium abundance in the older clusters. 
Loss of lithium is more severe the lower the mass of the star. 
Since the lithium depletion is driven by exposure to protons, 
even minor reductions of the 7 Li abundance below the value 
in zero age main sequence stars must be accompanied by 
very much larger reductions for 6 Li. 

Into this picture for lithium abundances of normal stars 
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Figure 10. Investigation of the effect of adopted Til abundance on the 6 Li/ 7 Li ratio for HD 82943. The dotted horizontal line is the 
observed Til abundance for HD 82943. The numbers are x 2 -values for the best fits for 1.88 eV case. 



must be woven the suggestion that main sequence stars 
may accrete one or more planets (or circumstellar material) 
which provide fresh lithium to the star's convective envelope 
with the largest effect being a replenishment of the 6 Li abun- 
dance. After accretion, the 6 Li/ 7 Li ratio in the star's atmo- 
sphere depends primarily on (i) the mass and composition 
of ingested material relative to the mass and composition 
of the stellar convective envelope, and (ii) the time elapsed 
since ingestion. Montalban & Rebolo (2002) have provided a 
useful exploration of these dependencies. Their calculations 
for main sequence stars consider the accretion of 1 to 10 M.j 
by main sequence stars of masses from about 0.7 Mq to 1.2 
Mq and of two compositions ([Fe/H] = 0.0 and 0.3). A ratio 
6 Li/ 7 Li ~ 0.10 immediately following accretion is achieved, 
for example, when a giant planet of mass 10 Mj is accreted 
by a 1.2 Mq main sequence. The post-accretion isotopic ra- 



tio decreases with decreasing mass of the star because the 
mass of a main sequence star's convective envelope increases 
to lower masses. Accretion increases also the 7 Li abundance, 
that is the 'total' lithium abundance. If lithium has not been 
depleted prior to accretion, the abundance increase is small 
(about 0.2 dex for accretion levels providing 6 Li/ 7 Li ~ 0.10), 
and, therefore, the total lithium abundance is not a sensitive 
indicator of the occurrence of accretion.^ 

Survival of 6 Li also determines its detectability. Within 
the constraints set by plausible assumptions about the ac- 

* Similarly, the beryllium abundance is a poor sensor of accre- 
tion. In fact, beryllium is largely undepleted in main sequence 
stars, even in stars showing lithium to be severely depleted. San- 
tos et al. (2002) measured Be abundances to find 'no clear differ- 
ence' between the abundances in stars with and without planets. 




Figure 11. Age estimation using the isochrones computed by Girardi et al (2000) for program stars with planets. 



cretion process and its benign effects on the structure of the 
star, the calculations show that accreted 6 Li (and 7 Li) is 
preserved for more than 3 Gyr in stars of mass greater than 
about 1.1 Mq. This limit corresponds to the effective tem- 
perature range T e g > 5900 K. Lithium-7 is predicted to be 
completely preserved, a prediction at odds with the lithium 
abundances for field and cluster stars (Figure 12). This dis- 
crepancy could be removed by using alternative models, see, 
for example, the predictions given by Proffitt & Michaud 
(1989). One can infer from observations summarized in Fig- 
ure 12 that accreted 6 Li will be depleted more quickly than 
predicted by Montalban & Rebolo's standard calculations: 
6 Li seems likely to survive for only about 1 Gyr in main se- 
quence stars with T eS > 6100 K. (For T eS ~ 6600 K, lithium 
is depleted to form the Boesgaard-dip [Boesgaard & Tripicco 
1986]), but both isotopes may be quite similarly depleted.) 

In the light of Figure 12 and Montalban & Rebolo's 
calculations, our failure to detect 6 Li in the young planet- 
hosting main sequence stars does not admit of a simple de- 
cisive conclusion about accretion of planets. Our measure- 
ments allow the following possibilities: (i) accretion has not 
occurred; (ii) accretion occurred but the accreted mass was 
insufficient to raise the star's 6 Li abundance above our detec- 
tion limit; and (iii) accretion occurred such that the post- 
accretion 6 Li abundance exceeded the detection limit but 
subsequent destruction of lithium depressed the 6 Li abun- 
dance again below the detection limit. Montalban & Re- 



bolo's predictions of the survivability of Li in stars with 
M > 1.1 Mq and our failure to see 6 Li may suggest that 
accretion can rarely add more than a Jovian mass of 'pri- 
mordial' material to the star. (Accretion of terrestrial rather 
than Jovian planets can effect an increase in the lithium 
abundances, but with a much smaller total mass for the ac- 
creted planets.) 

Effects of accretion of planets on surface compositions 
are not limited to the lithium abundance. Accretion of cir- 
cumstellar material has been invoked to account for the 
metallicity of planet-harbouring stars. There is now a gen- 
eral concensus that such stars are more metal-rich than sim- 
ilar stars without planets (Gonzalez et al. 2001; Santos et 
al. 2001; Reid 2002), but the origin of this difference is still 
a matter of debate. Gonzalez (1997) suggested two possi- 
ble explanations for the metallicity difference: accretion of 
H-poor matter (e.g., terrestrial planets, asteroids), and/or a 
bias in the discovery of planets in radial velocity surveys aris- 
ing from a metallicity- dependent migration of giant planets 
toward the star. Gonzalez (1998) added a third possibility: 
a high metallicity of the natal cloud as a prerequisite for for- 
mation of giant planets. Santos et al. (2001) favoured this 
latter possibility. Combinations of these possibilities may 
need to be considered; for example, if high metallicity clouds 
are favoured sites for stars with planets, such stars may ac- 
crete planets, and if these planets are terrestrial in nature, 
the metallicity of the star will be further increased. 
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Figure 12. The derived Li abundances for planet hosting stars are compared with the Li abundances of Pleiades, Hyades and M67 
clusters. Ages of the clusters are noted. The stellar masses of planet hosting stars are indicated by the size of the filled circles. 



Observational evidence in support of accretion as the 
origin of the high metallicity of planet-hosting stars is lim- 
ited. Detection of 6 Li in HD 82943 was claimed as supporting 
evidence. Our failure to confirm the presence of 6 Li weakens 
this claim, but, as we noted, accretion remains a possibility 
as long as 6 Li can be efficiently destroyed after the cessation 
of accretion. Accretion of giant planets will probably not 
change the surface composition of a star very much; giant 
planets are expected to have a composition similar to that 
of their host star, except for lithium and, perhaps, beryllium 
and boron. Since the composition of a terrestrial planet is 
far from a replica of the circumstellar material, the stellar 
composition will be changed by accretion of terrestrial plan- 
ets. One expects volatile elements (e.g., C, N, and O) to be 
reduced, and the abundance of non-volatiles (e.g., Fe-peak 
elements) to be increased by accretion. A signature of ac- 
cretion of terrestrial planets could be an abundance change 
correlated with an element's condensation temperature (T c ), 
the temperature at which an element condenses into solids 
as circumstellar gas cools in thermodynamic equilibrium. 
Gonzalez (1997) proposed this signature and searched for it 
in abundance analyses of two planet-hosting stars. Smith, 
Cunha, & Lazzaro (2001) noted that single (planet-less) 
metal-rich stars have a composition such that abundance 
differences, [X/H] (relative to solar values) for element (X) 
are quite well correlated with T c . These differences are, 



therefore, plausibly attributed to Galactic chemical evolu- 
tion, that is the net effect of nucleosynthesis by stars on the 
chemical composition of the interstellar medium. Stars host- 
ing planets show the same abundance - T c correlation at a 
given [Fe/H] indicating that accretion of terrestrial planets 
can have more than a minor role in setting the stars' sur- 
face composition (see Takeda et al. 2001). Smith et al. did 
draw attention to a sample of five planet-hosting metal-rich 
stars for which the [X/H] versus T c correlation appeared 
unusually strong and suggested that accretion of terrestrial 
planets might have occurred in these cases. Four of the five 
stars are in Table 1. 

Another example of a star which may have accreted 
dust, rocks, or terrestrial-like planets is HD 219542 A, a 
main sequence star with a distant main sequence compan- 
ion. Gratton et al.'s (2001) differential abundance analysis 
found HD 219542 A overabundant relative to HD 219542 B 
in non-volatile elements but not in the volatile elements. 
The former were overabundant by about 0.07 dex. Gratton 
et al. remark that accretion of a few earth masses of rocky 
material by HD 219542 A would suffice to account for the 
differential abundances, and might also explain the differ- 
ent lithium abundances in primary and secondary; lithium's 
condensation temperature is similar to that of silicon and 
iron. HD 219542 A is one of our programme stars and the Li 
abundance (Table 2) agrees well with that given by Gratton 



et al. (HD 219542 B is Li-poor, loge (Li) < 1.0, according 
to Gratton et al.) The lithium abundance of HD 219542 A 
is 'normal' for this young star. 

If accretion occured and added lithium, the 6 Li abun- 
dance corresponding to a 0.07 dex increase in the iron abun- 
dance would be approximately loge( 6 Li) = —0.1, if the ac- 
creted material contained a 'cosmic' abundance of lithium 
relative to iron. The corresponding 7 Li abundance is about 
log e( 7 Li) = 0.9 in the absence of lithium at the stellar sur- 
face prior to accretion. The observed abundance is consider- 
ably higher implying, as expected from Figure 12, that sub- 
stantial amounts of 7 Li were present when accretion occured. 
Combining the observed 7 Li abundance with the inferred 6 Li 
abundance from the mass of accreted rocky material gives 
6 Li/ 7 Li = 0.04. We measure 6 Li/ 7 Li = 0.03 ± 0.04 which 
is consistent with the inference but also with an absence of 
6 Li. Inspection of the observed profile gives a tantalising im- 
pression of a break in the red wing due to 6 Li (Figure 7). 
Higher S/N spectra will be sought. 



6 CONCLUDING REMARKS 

Our search for 6 Li in the atmosphere of planet-hosting star 
proved unsuccessful, despite the intriguing earlier report of 
substantial amounts of 6 Li in HD 82943, a star known to host 
two giant planets. A limit of approximately 6 Li/ 7 Li < 0.03 
was set for eight planet-hosting stars including HD 82943. 
Detection of 6 Li in a late F/early G-dwarf was recognized by 
Israelian et al. (2001) as a possible signature of accretion of 
planets by a star. Non-detection of 6 Li is less than watertight 
evidence against such accretion. Searches for planets should 
in the near future identify much larger numbers of planet- 
hosting stars. Accretion of circumstellar material may also 
occur onto stars without a planetary system. A continued 
search for 6 Li may yield a positive detection of 6 Li. 
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